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DARCY'S  LAW  FOR  FLOW  IN  POROUS  MEDIA 
AND  THE  TWO- SPACE  METHOD 


Joseph  B.  Keller 

Stanford  University 
Stanford,  California 


§1.  INTRODUCTION 

A  common  problem  in  engineering  and  science  is  to  derive 
simple  equations  governing  complicated  phenomena.  Often  com¬ 
plicated  governing  equations  are  known,  but  they  are  too  dif¬ 
ficult  to  analyze.  Therefore  it  is  desired  to  simplify  them 
in  such  a  way  that  the  phenomena  of  interest  will  still  be 
adequately  described,  while  finer  details  which  are  not  of 
interest  can  be  disregarded.  An  example  of  a  simplified  equa¬ 
tion  is  Darcy's  law,  which  describes  flow  of  a  viscous  fluid 
through  a  porous  medium.  The  more  complicated  equation  for 
the  3ame  phenomenon  is  the  Navier-Stokes  equation.  As  an  ex¬ 
ample  of  a  general  method  for  simplifying  equations,  we  shall 
show  how  to  derive  Darcy's  law  from  the  Navier-Stokes  equation 

Simplified  equations  are  often  called  "homogenized  equa¬ 
tions,"  and  the  procedure  of  replacing  the  original  equations 
by  them  is  often  called  "homogenization."  Traditionally, 
engineers  and  scientists  invent  homogenized  equations  from 
physical  considerations.  In  recent  years,  mathematicians  have 
developed  new  methods  for  deriving  them  from  the  original  com¬ 
plicated  equations.  These  methods  all  involve  averaging  of 
one  kind  or  another.  In  some  case  it  is  stochastic  averaging 
in  other  cases  it  is  spatial  averaging,  and  sometimes  it  is 
space-time  averaging. 
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The  introduction  of  averaging  depends  upon  the  fact  that 
the  problems  involve  two  widely  different  scales.  There  is  a 
long  scale  of  variation  which  is  of  interest  and  is  to  be  re¬ 
tained,  and  there  is  a  short  scale  of  variation  which  is  not 
of  interest  and  is  to  be  eliminated.  The  purpose  of  averaging 
is  to  eliminate  the  variations  on  the  short  scales.  Therefore 
the  region  over  which  a  spatial  average  is  taken  must  be  large 
compared  to  the  small  scale  and  small  compared  to  the  large 
scale. 

These  competing  requirements  have  traditionally  made  the 
derivation  of  simplified  equations  a  heuristic,  nonrigorous 
procedure.  However,  by  the  introduction  of  the  two-space 
method,  it  is  possible  to  convert  it  into  a  formal  and  rigor¬ 
ous  method.  In  fact,  whenever  apparently  contradictory  assump¬ 
tions  are  used  in  a  heuristic  analysis,  that  is  an  indication 
that  the  analysis  concerns  an  asymptotic  phenomenon.  Usually 
by  appropriately  recognizing  this  feature,  it  is  possible  to 
replace  the  heuristic  argument  by  a  systematic  and  formal 
asymptotic  analysis. 

The  mathematical  study  of  such  asymptotic  phenomena  has 
been  pursued  by  Spagnolo  [1],  DiCiorgi  and  Spagnolo  12], 
Sanchez-Palencia  [3,4],  Babuska  [5],  Bensoussan,  Lions,  and 
Papanicolaou  [6,7],  Larsen  [8,9],  McConnell  [10],  Keller  [11], 
and  others.  I  shall  employ  the  two-space  method  which  I  devel¬ 
oped  for  this  purpose  in  1973.  It  has  been  used  and  improved 
upon  by  Bensoussan,  Lions,  and  Papanicolaou  [6,7]  and  Larsen 
[8,9].  This  method  has  been  used  by  Burridge  and  Keller  [12] 
to  derive  the  equations  of  linear  poroelasticity  from  the  lin¬ 
earized  Navier-Stokes  equation  and  the  linearized  equations  of 
elasticity.  That  derivation  has  many  points  of  similarity  to 
the  present  one. 

Before  deriving  Darcy’s  law,  we  shall  consider  first  a 
simpler  example,  that  of  heat  conduction  in  one  dimension. 

The  thermal  conductivity  will  be  assumed  to  vary  rapidly  with 
position,  and  we  will  seek  a  less  rapidly  varying  effective 
conductivity.  In  Sec.  2  we  will  derive  it  from  the  explicit 
solution  of  the  heat  conduction  problem.  Then  in  Sec.  3  we 
will  rederive  it  by  the  two-space  method,  which  we  shall  later 
use  to  get  Darcy's  law. 
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§2.  EFFECTIVE  THERMAL  CONDUCTIVITY  IN  ONE  DIMENSION 

Let  us  consider  one-dimensional  heat  flow  along  a  rod  in 
which  the  thermal  conductivity  k  varies  very  rapidly  with  dis¬ 
tance  x  along  the  rod.  To  describe  rapid  variation  we  write 
k  =  k(x/e)  where  c  is  a  small  parameter.  Then  dk/dx  =  e_1k', 
so  dk/dx  is  large  when  e  is  small  even  though  k'  is  bounded. 

If  the  thermal  conductivity  also  varies  on  a  slow  scale,  we 
write  k  =  k(x,x/e).  Then  the  steady  state  temperature  distri¬ 
bution  along  the  rod  satisfies  the  equation 

4[k(x,x/e)  ^  u]  =  ^  9(x>  0  £  x  s  1  (2.1) 

Here  dg/dx  represents  a  heat  source,  which  has  been  written  as 
a  derivative  just  for  convenience,  and  we  also  assume  that 
g(l)  -  0.  As  boundary  conditions  we  will  suppose  that  one  end 
of  the  rod  is  kept  at  temperature  zero  and  the  other  end  is 
insulated,  so  we  have 

u  (o)  =  o  *  0  (2.2) 

The  solution  of  (2.1)  and  (2.2)  can  be  found  at  once  by 
integration.  By  using  the  assumption  that  g(l)  =  0,  we  get 
the  result 


u (x, e) 


(x') 

,x'/e) 


dx' 


(2.3) 


We  have  indicated  that  u  depends  upon  t,  which  is  a  measure  of 
the  small  scale  on  which  k  varies .  As  a  consequence  of  the 
small  scale  variation  of  k,  u  also  varies  on  this  small  scale. 
Our  goal  is  to  eliminate  this  small  scale  variation  of  u,  if 
possible,  by  considering  the  limit  of  u  as  e  tends  to  zero. 

We  shall  show  that  under  suitable  conditions  on  k,  u(x,e) 
has  a  limit  Uq(x), 


u„(x)  =  lim  u(x,e) 
u  e-i-0 


Furthermore,  uQ  is  given  explicitly  by 


uo(x)  58  j*  k^FT  dx' 


(2.4) 


(2.5) 
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Here  the  "effective  conductivity"  kQ(x)  is  defined  by 


Vx) 


=  lim 
£-*■0 

(■x+Ax 

dx' 

Ax 

Jx 

k(x,x'/e) 

=  lim 
e-»0 

e 

Ax 

t (x+Ax)/e  dy 

J  x/ £ 

(2.6) 


Our  hypothesis  about  k  is  that  the  limit  in  (2.6)  exists  uni¬ 
formly  in  x  and  is  independent  of  Ax. 

Before  proving  (2.5)  we  deduce  from  it  that  uQ(x)  satis¬ 
fies  the  equation 


s[‘oME'1W]-E’“  0<«S1  (2.7) 

and  the  boundary  conditions 
du  <D 

uQ(0)  «  0  — 33E~  *  0  (2.8) 


Thus  Uq(x)  satisfies  the  same  equations  as  does  u(x)  with 
k(x,x/e)  replaced  by  kQ(x).  This  justifies  the  name  "effec¬ 
tive  conductivity"  for  kQ(x).  We  note  that  kg(x)  is  the  har¬ 
monic  mean  of  k(x,x/e)  with  respect  to  its  second  argument. 

To  prove  (2.5)  we  consider  the  general  integral  I(x,e) 
defined  by 

I(x,e)  =>  [  f  (x*  ,x’/e)dx‘  (2.9) 

Jo 


We  let  Xj  *  jx/N,  j  *  0,  . ..,  N  -  1  where  N  >  1  is  an  integer, 
and  we  set  Ax  ■  x/N.  Then  we  rewrite  (2.9)  as  a  sum 


I (x,e) 


A 


X 

X 


j+1 

j 


flx’rX'/eiax' 


(2.10) 


Next  we  assume  that  f  is  continuously  differentiable  with  re¬ 
spect  to  its  first  argument  and  that  |fx|  ;  B  uniformly  in  c. 
Then  we  have 


|  f  (X*  ,X'/£)  -  f(XyX'/£)i 

X  <  x1  <  X  .  . 

xj  -  *  -  Aj+1 


|  (x’  -  Xj)fx(x. ,x’/e) !  5  BAX 

(2.11) 
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Here  x.  is  some  point  in  the  interval  x .  5  x .  s  x .  , , . 
3  333+1 

Now  we  can  use  (2.10)  and  (2.11)  to  write 


N-l  fx.  . 

I(x,e)  -  l  3+1  f (x.,x'/e)dx' 
j=l  Jx.  3 

3 

N-l  fX.  . 

s  l  3  \fix',x'/z)  -  f <x . ,x'/e) |dx‘ 

3=1  x •  3 

J  3 

<  NB(Ax)2  =  Bx2/N  (2.12) 


Upon  taking  the  limit  as  e  tends  to  zero  in  (2.12),  we  obtain 

N-l  r (x .+Ax)/e 

lim  I (x, e )  -  l  Ax  lim  t—  3  f(x.,y)dy 

e-0  j=l  e-0  flx  3x../e  3 

N-l  , 

=  lim  I (x, e)  -  l  Ax  f(x.)  s  Bx  /N  (2.13) 

e-*0  j=l  3 

Here  f (x)  is  defined  by  the  following  limit,  which  is  assumed 
to  exist  uniformly  in  x  and  to  be  independent  of  Ax: 

t (x+Ax)/ e 

fix)  =  lim  t —  f ( x , y ) dy  (2.14) 

e-0  Jx/£ 

Finally,  we  let  N  -»■  «  in  (2.13),  and  the  last  sum  tends  to  the 
_  2 

integral  of  f (x)  while  Bx  /N  tends  to  zero.  Thus  we  obtain 


lim  [X  f(x',xVe)dx' 
e-0  J0 


f 

I  n 


f ( x ' ) dx ’ 


(2.15) 


This  is  the  desired  result,  which  we  shall  summarize  as  a 
theorem. 


Theorem  Suppose  that  fx  is  continuous,  that  |fx(x,x/e)| 
s  B  uniformly  in  e,  and  that  the  limit  in  (2.14)  exists  uni¬ 
formly  in  x  and  is  independent  of  Ax.  Then  (2.15)  holds. 

The  result  (2.15)  might  be  called  an  integration  by  parts 
formula,  were  that  name  not  preempted,  since  part  of  the  in¬ 
tegration  is  done  first  in  (2.14)  and  the  other  part  later  in 
(2.15).  When  the  theorem  is  applied  to  (2.3)  it  yields  (2.5). 

From  the  method  of  derivation  of  the  theorem,  it  is  clear 
how  to  extend  it  to  the  case  of  more  than  one  rapidly  varying 
argument.  The  result  is 
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lim  [  f  (x' ,x'/e  ,x'A2, . . .  fX'/e^dx' 
e-*-0  J0 

=  f*  AlA?..-Anf(x*,y1,y2 . yn)dx'  (2.16) 

Here  A  ^  denotes  the  operation  of  averaging  with  respect  to 

the  variable  y..  This  extended  theorem  can  be  applied  to  (2.3) 
3  n 

in  which  k  =  k(x,x/£, _ ,x/e  ).  The  result  is  again  (2.5)  with 

kQ  given  by 

kQ (x)  *  Al" " "An  k (x.yjy  . . . ,yn>  (2.17) 

§3.  THE  TWO- SPACE  METHOD 

We  shall  now  solve  the  problem  of  the  preceding  section 
by  the  two-space  method,  which  is  applicable  to  more  general 
problems.  In  doing  so  we  shall  permit  the  source  term  to  be 
rapidly  varying  also.  Thus  we  replace  (2.1)  by 

3x(k(x,x/e) 3xJu  =  h(x,x/E)  0  <  x  s  1  (3.1) 


To  solve  this  equation  we  define  y  =  x/e,  write  u(x,e)  =  v(x,y,e) 
and  replace  3x  by  3x  +  e  *3  in  (3.1).  Then  (3.1)  becomes 

[e'23yk(x,y)3y  +  e'1(3xk3y  +  3yk3x)  +  3xk3x]v  =  h(x,y) 


(3.2) 


Now  we  seek  v  in  the  form 


2  2 

v(x,y,£)  =  Vg  (x,y)  +  ev^Xjy)  +  £  v2(x,y)  +  o(e  )  (3.3) 

We  substitute  (3.3)  into  (3.2)  and  equate  coefficients  of  the 
first  three  powers  of  e  to  obtain 


3yk3yv0  =  0 

3  k3  v,  =  -(3  k3  +  3  k3  )v_ 
yyl  xy  yx'O 


3yk3yv2  =  “<3xk3y  +  3yk3x)vl  “  3xk3xvO  +  h 


(3.4) 

(3.5) 

(3.6) 
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The  solution  of  (3.4)  with  yQ  arbitrary  is 

ry  -1 

vQ(x,y)  »  vQ(x,y0)  +  k (x,yQ) [3yvQ (x,yQ) ]  J  k  (x,y')dy' 

V°  (3.7) 

We  now  require  that  Vg(x,y)  be  a  bounded  function  of  y.  Then 
if  k  is  positive  and  bounded,  as  we  assume,  Vg  given  by  (3.7) 
will  not  be  bounded  unless  3yVg(x,yg)  =  0.  Since  y^  is  arbi¬ 
trary  this  implies  that  Vg  is  independent  of  y: 

Vg  =  Vg(x)  (3.8) 

Upon  using  (3.8)  in  (3.5)  and  solving  for  v^,  we  get 

v1(x,y)  =  v1(x,yQ)  -  (/  -  Y0)3xvo*x*  +  k(x'Vo^3yvl  ^x,y0^ 


3xv0(x)l  jY  k-1 (x,y 1 )dy’ 


We  next  impose  the  requirement  that  v^  be  a  bounded  func¬ 
tion  of  y.  Then  upon  dividing  (3.9)  by  y  -  yQ  and  letting 
y  -*  oo,  we  get 


3xv0(x)  =  k(x'y0)t3yVl(x'y0>  +  3xv0(x)lk0  (x) 


(3.10) 


Here  kQ  (x)  is  defined  by  the  following  limit,  which  we  assume 


to  exist  and  to  be  independent  of  yQ: 


k”1(x)  =  lim 

y-*-oo 


y  -  y0 


i—  fY  k-1  (x,y  1 ! 

y0  JV. 


(3.11) 


Now  we  solve  (3.10)  for  k3yv^  and  substitute  the  result 
into  (3.6),  which  we  rewrite  as 

3yk3yv2  +  3yk3xvx  =  -3xtkg (x) 3xvQ (x) ]  +  h(x,y)  (3.12) 

integration  of  (3.12)  from  yQ  to  y  yields 

(k3  v2  +  k3xVl)  Y  =  -(y  -  y0)3x[kg(x)3xv0(x)] 


[Y  h  ( x ,  y  ’ )  dy  ’ 

;yo 


(3.13) 
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Dividing  (3.13)  by  y  -  y^  and  letting  y  °°  yields,  in  view  of 
the  assumed  boundedness  of  terms  on  the  left  side  of  (3.13), 

3x[kO(x)3xvO(x)1  =  S(x>  (3.14) 

Here  h(x)  is  the  average  of  h(x,y)  with  respect  to  y. 

The  result  (3.14)  is  the  desired  equation  for  vQ(x),  the 
leading  term  in  u(x,e)  =  v(x,y,e).  The  effective  conductivity 
kg(x)  given  by  (3.11)  is  the  same  as  that  defined  by  (2.6)  in 
Sec .  2 . 

§4.  FLOW  IN  POROUS  MEDIA 

Let  us  now  consider  the  flow  of  a  compressible  viscous 
fluid  through  a  rigid  porous  medium.  The  pore  "diameter"  is 
assumed  to  be  small  compared  to  the  macroscopic  scale  of  the 
medium.  We  shall  denote  by  e  the  ratio  of  the  pore  diameter 
to  the  macroscopic  scale,  and  by  D£  the  interior  of  all  the 
pores.  Then  the  equations  governing  the  fluid  flow,  which  are 
the  Navier-Stokes  equation,  the  continuity  equation,  and  the 
equation  of  state,  hold  in  De .  These  equations  are 


2  1 

p(3fcu  +  u*Vu)  =  -Vp  +  p  (V  +  Vtf’)u  +  f 

(4.1) 

3fcp  +  V • (pu)  =  0 

(4.2) 

P  =  P(p) 

(4.3) 

Here  p,  p,  u,  and  p  are  the  fluid  density,  pressure,  velocity, 
and  viscosity  coefficient,  respectively,  while  f  is  the  exter¬ 
nal  body  force  per  unit  volume.  In  addition  to  these  equations, 

we  have  u  =  0  on  3D  . 

e 

Our  goal  is  to  derive  simplified  equations  for  the  fluid 
motion  by  taking  advantage  of  the  smallness  of  e.  To  this 
end  we  first  introduce  y  and  t  defined  by 

y  =  x/e  t  =  t /e3^2  (4.4) 

Then  we  write  u,  p,  p,  p,  and  f  in  the  forms 
1/2 

u  =  e  u(x,y,t,e)  p  =  p(x,y,x,e)  p  =  p(x,y,t,e) 

f  =  f(x,y,t,e)  p  =  e3//2p  (4.5) 


(4.16) 
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PQ  =  p  0 (x  #  t )  =  p[p0<x,T)J  (4.17) 

Next  from  (4.1)  we  find  that  V^'u q  =  -Pg^a^Pg,  and  from  (4.17) 
we  conclude  that  V^.u^  is  independent  of  y.  Therefore, 

VyVy*u0  =  0,  so  this  term  can  be  omitted  from  (4.13). 

Now  we  consider  (4.11)  and  (4.13)  as  a  system  of  equations 
for  Uq  and  p^  as  functions  of  y  and  i ,  in  which  x  occurs  as  a 
parameter.  In  these  equations  we  view  pq(x,t)  and  fg(x,x)  - 
Vxpq(x,t)  as  a  given  time  dependent  density  and  a  given  time 

dependent  body  force,  respectively.  Both  of  them  are  indepen¬ 

dent  of  y.  We  seek  a  solution  Ug,  p^  of  these  equations  with 
the  normal  component  of  Ug  vanisning  on  the  boundary  of  the 
domain,  which  may  depend  upon  x.  If  the  medium  is  spatially 
periodic  in  y,  we  require  the  solution  to  have  the  same 
periodicity;  otherwise,  we  require  it  to  be  a  bounded  function 
of  y.  In  either  case  we  assume  that  the  solution  is  unique. 
Then  we  can  write  it  as  a  functional  of  p.  and  f.  -  V  p.  in 
the  form 

u0(x«y/i)  =  U[x»Y.T»P0.f0  "  VxpQ]  (4.18) 

P1(x,y,x)  =  P [x,y , i , pQ , f 0  -  VxpQ]  (4.19) 

The  functionals  U  and  P  involve  their  argument  functions  only 
for  t ’  ST. 

We  now  define  the  average  of  any  function  of  y  by  inte¬ 
grating  it  with  respect  to  y  over  a  large  domain  D  of  the 
fluid,  dividing  the  integral  by  the  volume  V  of  D,  and  letting 
D  and  V  tend  to  infinity.  In  particular,  by  averaging  (4.18) 
and  (4.19)  we  obtain 

uq(x,t)  =  U[x,x,p0,f0  -  VxpQ] 

=  lim  i  [  U[x,y,T,p0,fQ  -  VxpQ]dy  (4.20) 

V-*°o  *  D 

Pi^'T)  =  PU,T,P0,f0  -  Vxp0) 

=  lim  i  P  [x,y,T  ,p-,f-  -  Vxp,Jdy 
V+oo  *  D 


(4.21) 
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These  two  relations  express  uQ  and  as  functionals  of  VxPg 
and  Pq,  so  together  they  represent  a  generalization  of  Darcy's 
law  to  nonlinear  time  dependent  compressible  flows.  We  shall 
examine  them  further  in  the  next  section. 

§  5 .  DARCY ' S  LAW 

Relations  (4.20)  and  (4.21)  can  be  used  to  derive  an 
equation  for  pg(x,x).  To  obtain  it  we  average  (4.14)  in  the 
manner  just  described.  Integration  of  the  term  V  • (pgU^  + 
p^Ug)  yields  a  surface  integral  over  the  "ends"  of  the  pores, 
since  uQ  and  u^  vanish  on  the  rigid  boundaries.  When  divided 
by  V  this  integral  tends  to  zero,  and  then  (4.14)  leads  to 

3Ttpp(p0)Pi]  +  vx‘tPo^O]  =  0  (5*1) 

In  (5.1)  we  have  used  (4.15)  to  replace  p^  by  PpP^*  Now  we 
use  (4.20)  and  (4.21)  for  Ug  and  p^  in  (5.1)  to  obtain 

3T{Pp(Po)P^,T,p0,fo  -  VxPg]} 

+  Vx*{pgU[X,T,Pg,fg  -  VxPg] }  =  0  (5.2) 

Equations  (5.2)  and  (4.17)  are  a  pair  of  equations  for  pQ(x,t) 
and  pg(x,t).  If  (4.17)  is  used  to  eliminate  pQ,  the  result  is 
a  single  equation  for  Pg(x,x). 

We  shall  now  consider  some  special  cases  of  (4.20)  and 
(4.21)  which  are  closer  to  the  usual  form  of  Darcy's  law. 

First,  for  an  incompressible  fluid  of  constant  density  Pg,  we 
need  not  indicate  the  argument  pQ  in  (4.20),  which  then  becomes 

Ug(X,T)  =  U[X,T,fg  -  VxPg]  (5.3) 

Furthermore,  p^  =  0  so  p^lpgJP  =  0,  and  then  (5.2)  simplifies 
to  the  following  equation  for  pg(x,x): 

Vx*U[x,x,fg  -  VxPg]  =  0  (5.4) 

In  the  steady  incompressible  case  U  is  just  a  function  of 
fg  -  7xPg,  rather  than  a  functional,  and  it  is  independent  of 
x.  Then  (5.3)  and  (5.4)  simplify  accordingly,  and  (5.3)  is  a 
nonlinear  form  of  Darcy's  law  for  steady  incompressible  flow. 
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In  the  case  of  steady  compressible  flow,  (4.20)  and  (4.21) 
become 

uQ(x)  =  U[x,pg(x) ,fQ(x)  -  Vxp0(x)]  (5.5) 

(x)  =  P[x,pQ(x) , f 0 ( x)  -  VxP0(x)]  (5.6) 

Here  U  and  P  are  ordinary  functions,  and  not  functionals,  uo 
(5.5)  is  a  nonlinear  Darcy  law  for  steady  compressible  flow. 

The  corresponding  equation  for  pQ  is,  from  (5.2), 

Vx‘ {pO0[x'°O'fO  "  Vxp0])  =  0  (5.7) 

Finally,  we  shall  consider  the  case  in  which  pq(x,t)  is 
close  to  a  constant  density  m  and  fg  -  VxpQ  is  close  to  zero. 
Then  we  can  linearize  the  functionals  U  and  P.  In  doing  so, 
it  is  helpful  to  return  to  Eqs .  (4.11)  and  (4.13)  by  means  of 
which  U  and  P  are  defined.  The  linearized  forms  of  these 
equations  are 

3tP0  +  mVy»uQ  =  0  (5.8) 

m3Tu0  =  -Vypx  +  u(V^  +  i  VyVy*)uQ  +  fQ  -  VxpQ  (5.9) 

We  seek  the  bounded  solution  uQ,  p^  of  these  equations 
which  has  vanishing  velocity  on  the  rigid  boundaries  of  the 
pores.  This  solution  is  a  linear  functional  of  the  inhomoge¬ 
neous  terms  3tpg  and  fQ  -  7xPg.  Therefore  we  write  it  in  the 
form 

uq (x,y,r)  =  A (x,y , t ) (f g  -  ^XPQ)  +  B(x,y,T)3TP0  (5.10) 

Pj_(x,y,T)  =  C  (x,y ,  t  )  (f  g  -  VxpQ)  +  E (x,y , t ) 3tp0  (5.11) 

Then 

Ug(X,T)  =  A  (x,  T  )  (f  g  -  VxPg)  +  B  (  X  ,  T  )  3  t  P  g  (5.12) 

Pj^fx,!)  =  C(x,x)(fg  -  VxPg)  +  E(x,t)3tp0  (5.13) 

These  results  are  the  linearized  forms  of  the  generalized 
Darcy  law  for  time  dependent  compressible  flows. 
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The  linearized  form  of  (5.1),  in  view  of  (5.12)  and  (5.13), 
is 

Pp(*)3T[C(f0  -  VxpQ)  +  E3tp0] 

+  mVx-[A(f0  -  VxpQ)  +  B3tP01  =  0  (5.14) 

Here  it  is  the  constant  unperturbed  pressure  corresponding  to 
the  density  m. 

In  the  steady  case  (5.12)  simplifies  to  the  usual  form  of 
Darcy  1 s  law 

Ug  (x)  =  A(x)[fg(X)  -  VxPg(X)]  (5.15) 

Here  A(x)  is  just  a  matrix  which  transforms  as  a  tensor.  The 
steady  form  of  (5.14)  is  a  linear  equation  for  pQ (x) : 

Vx-{A(x)[fQ  -  VxpQ]}  =  0  (5.16) 

For  an  incompressible  fluid  (5.12)  simplifies  to 
UQ(X,T)  =  A  (X ,  T  )  ( f  Q  -  VjjPq)  (5.17) 

Since  pp  =  0  in  this  case,  (5.14)  becomes 

V{S(fo  -  Vo)}  "  0  (5-18) 

In  the  steady  incompressible  case,  (5.17)  reduces  to  the  usual 
Darcy  law  (5.15).  When  the  medium  is  macroscopically  isotropic, 
A (x)  is  a  scalar. 


§6.  CONCLUSION 

After  illustrating,  on  one-dimensional  heat  conduction, 
the  two-space  method  for  deriving  simplified  equations,  we 
have  applied  it  to  flow  in  porous  media.  Then  we  were  able 
to  derive  Darcy's  law  and  various  generalizations  of  it.  To 
determine  the  coefficients  in  this  law  and  its  generalizations, 
it  is  necessary  to  solve  a  difficult  problem  which  we  have  not 
considered.  That  problem  involves  the  detailed  pore  configura¬ 
tion.  Our  goal  was  to  determine  the  form  of  the  law  and  to 
characterize  the  coefficients  in  terms  of  the  solution  of  a 
specific  problem.  In  the  heat  conduction  problem  we  could 
actually  find  the  coefficient  because  the  problem  was  so  simple. 
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